Detection and characterization of multipartite entanglement in optical lattices 
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We investigate the detection and characterization of entanglement based on the quantum network 
introduced in [Phys. Rev. Lett. 93, 110501 (2004)] for different experimental scenarios. We first 
give a detailed discussion of the ideal scheme where no errors are present and full spatial resolution 
is available. Then we analyze the implementation of the network in an optical lattice. We find 
that even without any spatial resolution entanglement can be detected and characterized in various 
kinds of states including cluster states and macroscopic superposition states. We also study the 
effects of detection errors and imperfect dynamics on the detection network. For our scheme to be 
practical these errors have to be on the order of one over the number of investigated lattice sites. 
Finally, we consider the case of limited spatial resolution and conclude that significant improvement 
in entanglement detection and characterization compared to having no spatial resolution is only 
possible if single lattice sites can be resolved. 

PACS numbers: 03.75.Gg, 03.67.Mn 



I. INTRODUCTION 

Multipartite entanglement is an essential ingredient for 
complex quantum information processing (QIP) tasks, 
such as quantum error-correction [lj, multi-party quan- 
tum communication and universal quantum compu- 
tation in the one-way model 0- It has been generated 
between photons and in a controlled way both in ion 
traps |5| and in experiments with Mott insulating states 
of neutral bosonic atoms in optical lattices |(| ■ For 
photon and ion trap experiments full state tomography 
|j, |S| impressively proved the existence of multipartite 
entanglement with few particles. However, its unambigu- 
ous detection in optical lattice setups has so far not been 
possible. The measurements implemented in the actual 
experiments only provided information on the average 
density operator of each individual atom (representing a 
qubit), but not on the density operator of the compos- 
ite system |6{. Full state tomography in optical lattices 
seems to be infeasible due to the large number of involved 
atoms and thus other methods are required to prove the 
existence of entanglement. Recently, a quantum network 
that can detect multipartite entanglement in bosons, and 
can be realized in optical lattices, was proposed 0]. 

This quantum network is made up of pairwise beam- 
splitters (BS) acting on two identical copies of a multi- 
partite entangled state p\2... n of n atoms, and allows the 
determination of the purity of pn... n and of all its reduced 
subsystems. The purities provide us with a separability 
criterion based on entropic inequalities that, if vio- 
lated, indicate that the state pn...n is entangled. This 
nonlinear entanglement test can be implemented experi- 
mentally with a number of copies of pii... n independent 
of n, rather than the exponentially large number of copies 
needed for a full state tomography, and is more powerful 
than the other common experimental methods, such as 
Bell inequalities [lfj and entanglement witnesses . In 
fact, the nonlinear inequalities are known to be strictly 
stronger than the Bell-CHSH inequalities [l^. Further- 



more, if it may be assumed that the entangled state is 
characterized by a few parameters only it is often pos- 
sible to determine them from the results obtained from 
the entanglement detection network. 

The main aim of this paper is to analyze the opera- 
tion of the entanglement detection network under realis- 
tic experimental conditions. After introducing the ideal 
scheme we consider the particular realization of the net- 
work in an optical lattice. We discuss three main sources 
of imperfections: limited or no spatial resolution, errors 
in the BS operations occurring with a probability q, and 
failure to detect an atom with probability p. In current 
experiments none of these errors can be suppressed com- 
pletely and thus it is important to find out whether their 
presence still allows the unambiguous detection and char- 
acterization of multipartite entangled states. Our main 
results will be that (i) even without spatial resolution it 
is possible to unambiguously detect and characterize en- 
tanglement in multipartite states like cluster states and 
macroscopic superposition states, (ii) the effects of BS 
and detection errors can be made arbitrarily small by in- 
creasing the number of runs of the entanglement detec- 
tion network, and (iii) the required number of experimen- 
tal runs for achieving this is reasonably small if kq < 1 
and kp < 1, where k < n is the number of particles in 
the subsystem whose reduced purity is being measured. 
Finally we will find that (iv) single site resolution is nec- 
essary to achieve significant improvement compared to 
having no spatial resolution at all. 

The paper is organized as follows: in Sec. ^ we in- 
troduce the entropic inequalities and discuss how multi- 
partite entanglement can be detected and characterized 
by them. Then we describe the entanglement detection 
network and show how it can be realized in optical lat- 
tices and how the main sources of imperfections arise. In 
Sec. IIIII we investigate the detection and characterization 
of multipartite entangled states if no spatial resolution 
is available, and in Sec. IIVI we analyze the effect of the 
dominant experimental errors. We also discuss the case 
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of limited spatial resolution in the measurements. Finally 
we summarize our results in Sec. ]V\ 



II. ENTANGLEMENT DETECTION AND 
CHARACTERIZATION IN OPTICAL LATTICES 

In this section we first present the entropic inequalities 
which are used to detect multipartite entanglement. We 
show how they can be utilized to characterize multipar- 
tite entanglement if the entangled state is described by 
few parameters only. Then we introduce our entangle- 
ment detection network and discuss its implementation 
in optical lattices. Finally we study the main types of 
imperfections and resulting errors that we expect to be 
present in this experimental realization. 



A. Entropic inequalities 

Our network uses the entropic inequalities introduced 
in for multipartite entanglement detection. These in- 
equalities provide a set of necessary conditions for sepa- 
rability in multipartite states. Consider a state pi23...n 
of n subsystems. If pns...n is separable, then 



Pl23.. 



(i) 



where p^ is a state of subsystem j, and Cn = 1. The 
purity Tr {p 2 B } of the reduced density operator p B with 
B C {1,2, ... ,n} is less than or equal to the purity of 
any of its reduced density operators: 



Tr {p 2 A } > Tr {p 2 B } for all A<ZB. 



(2) 



We can use these inequalities to give a relation between 
the average purities of all reduced density operators of a 
given number k of subsystems defined by 



Tr {^)} 



E Tr K) 



(3) 



\B\ = k 



where B is summed over all combinations of k subsys- 
tems. For completeness we define Tr |p( )} = ^' F rom 
Eq. J2J we find that for any separable state 



Tr 



{p\ k) }>^{p\ kl) }ior allfc<fc'. 



(4) 



We will now consider how Eqs. J5J can be used to 
detect entanglement and to characterize entangled states. 



1. Entanglement detection 

Any state p that violates any of the inequalities 
Eqs. J21EJ) is entangled. If pi2d,... n is separable and pure, 



Tr { ( 0i23...„} = 1 and all the above inequalities become 
equalities; since a state with all Tr {p%} = 1 is nec- 
essarily a product state, all pure entangled states vio- 
late Eqs. J21 EJ- All pure entangled states can hence 
be detected by comparing Tr {p 2 23 «} with any other 

Tr {p%} or Tr |p(a,)| unlike entanglement witnesses or 

Bell inequalities where different states often require dif- 
ferent witnesses. From the Schmidt decomposition of any 
pure state into two disjunct subsystems A and B with 
A U B = {1,2... n} we find Tr {p 2 A } = Tr {p 2 ,} and 

hence Tr = Tr |p| n _ fe) | for all pure p. 

Mixed entangled states do not always violate Eqs. J5J 
0}, but since Tr {p%} is continuous in p a sufficiently 
small amount of noise added to a pure entangled state 
will leave it still violating the inequalities. In the ex- 
amples studied in this paper the noise level at which 
the inequalities no longer detect entanglement is a large 
fraction of the level at which entanglement ceases to be 
present. The permissible range values for the purities of 
any state p is 



2- fc <Tr{p|}<l or 2" fc < Tr {pf k) } < 1, (5) 

where k = \B\. The minimum is attained by the maxi- 
mally mixed state p = 2~ n I (with I the identity operator) 
and the maximum by any pure product state. 

Our test can often verify that a state is truly n-particle 
entangled and not just a collection of small entangled 
subsystems, because Tr {p%} < Tr {p 2 } always indi- 
cates that B is entangled with the rest of the system 
and cannot be caused by entanglement within B itself. 
Furthermore our test is insensitive to local unitary op- 
erations, which cannot alter the entanglement structure, 
since they do not change the purities. It is, however, 
in most cases sensitive to local measurements destroying 
entanglement as these generally do change the purities. 
If the tested state p is of a known form and characterized 
by few parameters only the kind and degree of violation 
of Eqs. © can be used to determine these parameters 
and thus characterize the state. 



2. Characterization of entanglement: Macroscopic 
superposition states 

Macroscopic superposition states | 7„) defined by 



l7n> 



|or + (7|o) + VT~M T |i))" 

+ 7" + T 



(6) 



with a single complex parameter 7 satisfying I7I < 1 
arise naturally in several systems such as BECs [l3| and 
superconductors Q. The quantity n(l — I7I 2 ) has been 
suggested as a measure of the effective size of the state, 
as in some respects | 7„) is equivalent to a maximally 
entangled state of n(l — I7I 2 ) particles The purities 
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of 1 7„) are given by 



Tr {pi} = 



2 + 27 fe 7 fe + 27™- fc 7 n - fe 



(2 + 7™ + 7™) 2 
47™ + 47™ + 7 2 " + 7 
(2 + 7™ + 7™) 2 



-2n 



(7) 



where k = \B\. These purities violate Eq. (J2J), showing 
that I 7„) is entangled, for all I7I < 1 and measuring the 
purities allows determining 7. For 7 = the state | 7^) is 
a maximally entangled state which is pure (Tr {p 2 } = 1) 
while all its subsets have purity Tr {p%} = 1/2. Up to 
local unitaries it is the only state with these purities, so 
can be unequivocally identified by our test. 

We demonstrate the effect of noise on our method by 
studying individual dephasing of each qubit. Such de- 
phasing noise maps a state according to 



£ - 



n 



(8) 

where cr ZJ is a phase flip applied to qubit j, A is the set 
of qubits which undergo a phase flip and < d < 1 is 
the decoherence parameter. Applying the map Eq. (JHJ) 
to I (7 = 0)„) we find for the resulting density operator 

p = |(ior<or+]i>"<i n 



(i-dy 



(|o)"(ir + |i)"(oD, (9) 



and thus Tr {p 2 } = (1 + (1 - d) 2n )/2 while Tr {p 2 B } = 
1/2 for all subsets B. Therefore entanglement is detected 
whenever it is present. This is not generally the case as 
can easily be seen by looking at a Werner state p — (1 — 
d) I (7 = 0)„) ((7 = 0)„ 1 + 2— "d I for which entanglement 
is detected if d < 1 — (2™ _1 + l) -1 / 2 . However, e^s. in the 
case n = 2 this state is entangled iff d < 2/3 [l6|, while 
our test works only for d < 1 — 1/V3- 



3. Characterization of entanglement: Cluster states 

Cluster and graph states are multipartite entangled 
states which form the basic building blocks of the one- 
way quantum computer 0|. We consider cluster- like 
states I 4> n ) defined by 



1 



11(1°), 



1 



2™-l 



/on Z — ^ 



(10) 

where c(x) is the number of occurrences of the sequence 
01 in the n-bit binary number x. These states have al- 
ready been created in optical lattices and represent 
a cluster state for <f> = it. Current methods for detect- 
ing them essentially perform a tomographic measurement 



of the average single particle density matrix, which goes 
from pure (| 0) + 1 1))((0 | + (1 1)/2 at <j> = to maximally 
mixed 1/2 at <fi = ir and back again at <j) = 27r. This 

method thus yields one measurement Tr |p^i)| relating 
to entanglement and two measurements relating to local 
unitaries. 

Our network permits the measurement of further cor- 
relations and it does not need the assumption that all 
atoms have the same single particle density matrix. It 
therefore allows a better characterization of \4> n )- For 
any <fi which is not an integer multiple of 2tt 7 \4> n ) is a 
pure state with no separable subsystems, and hence for 
any subset B we have Tr {p B } < 1 = Tr {p 2 } as can be 
seen from Fig.Q] For creating the states | 4> n ) each qubit 
only needs to interact with their two nearest neighbors 
except for the two extremal atoms 1 and n which will 
interact with only their one neighbor. Because of this 
the reduced purity of a subset B is determined by the 
boundary between it and the rest of the system. Subsets 
of different sizes but with the same boundary structure 
have the same purity. Several examples of these purities 
are shown in Fig. E]as a function of <j>. For example, all 
sets of two or more adjacent atoms, located anywhere in 
the row that do not include either extremal atom have 
the same purity (1 + cos 2 (</)/2)) 2 /4 (dash dotted curve in 
Fig. ^) which is independent of n. The degree of viola- 
tion of Eq. J2Jl varies smoothly with <j) and measuring the 
various different purities allows to determine 4> (up to its 
sign). We will now introduce a quantum network which 
detects violations of Eq. J2J and later, in Sec. 1 1 1 1 1 show 
how violations of Eq. |0J can be detected even without 
achieving spatial resolution of the different subsystems. 



B. Multipartite Entanglement Detection Network 

A family of quantum interferometric networks that 
directly estimates Ti{p 3 }, s = 2,3,4,... for any p 
from s copies of p was introduced in |17|. These net- 
works rely on the controlled-shift operation C — be- 
tween the different copies, where \ 4>\) \ $2) ■ • ■ | 4> s ) = 
I <j) s ) I <pi) . . . I (j) s -i) for any | fa) with i = 1, 2 . . . s. For 
the particular case of s = 2 the value of Tr {p 2 } is di- 
rectly related to the probability of projecting p (g) p into 
its symmetric and antisymmetric subspaces. The val- 
ues of the 2™ different purities associated with pi23... n 
can hence be determined from the expectation value of 
the symmetric and antisymmetric projectors, on each dif- 
ferent pair of reduced states ps = Tr a(pi23. ..«.), where 
, A U B = {l,2,...n}. These expectation values can be 
measured without resorting to the implementation of the 
three-qubit C-Swap gate C — if two identically pre- 
pared ID rows of n qubits which are represented by 
bosonic particles are coupled via pairwise BS, as shown 
in Fig.H0]. 

The BS in the j-th column projects the symmetric (an- 
tisymmetric) part of the density operator pj onto doubly 
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FIG. 1: Effect of varying <f> on Tr {p|} where B is any one 
atom not at an end (dotted), any two atoms not at ends and 
with at feast two others between them (dashed), any two or 
more consecutive atoms not inciuding an end (dash-dotted), 
any one or more consecutive atoms inciuding one end (soiid). 
These purities are independent of n. For n > 6 correspond- 
ing subsets are e.g. B = {2}, {2, 5}, {2, 3, 4, 5}, {1, 2, 3, 4, 5}, 
respectiveiy. Hence the state is entangied if the purities as 
fisted above are not in decreasing order that is for aff 4> ex- 
cept integer multipies of 2n. 
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FIG. 2: Entangfement detection network acting on two rows 
/ and II of identicai qubits. The two rows of n qubits are 
prepared in the (possibiy entangfed) state pi23...n <8> pi23...n 
and performing the pairwise BS foflowed by measuring the 
iattice site occupancies reaiizes the entangfement detection 
network by projecting each operator pj into its symmetric 
and antisymmetric subspaces. 



Tr {pj} = P+++ + P ++ _ + P+_+ - 

-P-4-+ - P-4- - P—A 



P+— (16) 
- P 



C. Realization in Optical lattices 



(singly) occupied sites 0] (for details see Appendix El . 
Two qubits in column j and in state pj ® pj will thus 
end up in the same site (+) or in different sites (-) with 
probabilities 

P® = ~Tr {(I ± V^) Pj ® Pj }=\± \lr {p]} . (11) 

Here S± = (I±V < - 2 ))/2 is the symmetric / antisymmetric 
projector. By distinguishing doubly occupied sites from 
singly occupied ones we can thus determine the purity of 
Pj- 

Extending this two-qubit scenario to the general case 
of two copies of a state of n qubits undergoing pairwise 
BS (see Fig. EJ) we obtain @ 

P±i± 2 -±„ = Tr |JJ y Pi2...n<8pl2...n| ■ (12) 

Inverting the linear equation Eq. (|12|l the impurity of 
any subset of atoms B is given by twice the probability 
of having an odd number j B of antisymmetric projections 
in subset B 

Tr {pi} = P{j B even) - P(j B odd) = 1 - 2P(J B odd). 
For example, for n = 3: 

Tr {p 2 123 } = P+++ + P+— + P-+- + P—+ (14) 

~P _ P+-+ ~ P++- _ P-++- 

Tr {p 2 12 } = P+++ + P++- + P—+ + P— (15) 

— P-++ - P-+- - P+-+ - Pa > 



We consider one sheet in the x — y plane of an ultracold 
two-species bosonic gas confined in a 3D optical lattice 
sufficiently deep so that the system is in a Mott insulat- 
ing state with exactly one atom per lattice site [19( ■ Two 
long lived internal states of each atom | a) and | b) repre- 
sent a qubit, i.e. for a site with row index I and column 
index j we define two basis states (ap )t|vac) = | 0); and 

(6p' l )t|vac) = | 1);. Here |vac) is the vacuum state and a 
(b) is the bosonic destruction operator for an atom in in- 
ternal state | a) (| b}). We assume that starting from this 
Mott state identical multipartite entangled states p\2z... n 
are created in each row while different rows remain un- 
corrected. This can e.g. be achieved by state selective 
cold controlled collisions between atoms in neighboring 
columns 0, |2(| • The above entanglement detection net- 
work can be realized in this setup to study the entangle- 
ment properties of pva...n- We will first discuss the im- 
plementation of the BS by coupling pairs of rows along 
the x-direction and then investigate methods to distin- 
guish doubly from singly occupied sites. Note that the 
intrinsic parallelism of the 3D optical lattice will allow to 
run many copies of these networks in the lattice at once. 
By exploiting this parallelism one can obtain an estimate 
of the desired projection probabilities and test the viola- 
tion of Eqs. (0 in a single or few experimental runs 
only. 

1. The pair-wise Beam Splitters 

The dynamics of the atoms in the optical lattice is gov- 
erned by the Bose-Hubbard Hamiltonian (BHM) p^.l20|. 
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We assume that any hopping along the y and z direc- 
tions is suppressed by a sufficient depth of the lattice in 
these directions. The hopping in x-direction is controlled 
by dynamically varying the corresponding lattice depth 
Vox- Since Vq x is proportional to the laser intensity it 
can easily be changed in the experiment. The pairwise 
BS requires that rows be coupled pairwise (this can be 
achieved by using a superlattice of twice the period) and 
thus we only need to consider one such pair labelling the 
two rows by I and 77, respectively, as shown in Fig. 

The Hamiltonian describing the dynamics of the atoms 
in these two rows can be written as a sum Hbhm — 
Hbs + Hjj where Hbs is due to hopping along the x 
direction and Hjj is due to the repulsion between two 
atoms occupying the same lattice site 0, |2£j . The two 
contributions are given by 

n 

Hbs = -J^f^n+b^b^ +h.c), (17) 

1=1 JI 3=1 

where J is the hopping matrix element and U is the onsite 
interaction energy. The parameters U and J depend on 
the parameters of the trapping lasers, and their ratio can 
be varied over a wide range, on time scales much smaller 
than the decoherence time of the system, by dynamically 
changing the depth Vb of the optical lattice [ljj, |2(j ■ 

The BS dynamics is perfectly realized in the non- 
interacting limit U = by applying Hbs for a time 
tbs = 7r/(4J) (for details see Appendix 0. However, in 
practice it is impossible either to control J perfectly ac- 
curately or to completely turn off the interaction U, and 
these imperfections cause the symmetric component to 
have a non-zero probability qb s of failing to bunch which 
is given by 

TT 2 [5J\ 2 1 [U\ 2 

^ir(-) + -6{i) • (18) 

Here 5J describe the fluctuations in J (for details see 
Appendix 1X1) . If the fluctuations 8 J occur from run to 
run rather than from site to site qbs can be interpreted 
as a statistical random variable. We will discuss how to 
correct this error in Sec. IIVI 

2. Measuring the lattice site occupation 

Recently, a method that uses atom-atom interactions 
to distinguish between singly and doubly occupied sites 
was demonstrated experimentally pi) . However, a sim- 
plified alternative method where rapid same-site two- 
atom loss is induced via a Feshbach resonance and the re- 
maining singly-occupied sites are detected suffices for our 



purpose. The detection of the remaining singly-occupied 
sites is achieved by measuring the atomic density pro- 
file after the atoms are released from the lattice. A 
single Feshbach resonance will cause the loss of either 
\aa), \ab) + |6a), or | bb). Hence, in order to empty 
doubly-occupied sites in all three states we can either 
turn on consecutively three separate Feshbach resonances 
or change the internal state of the pairs of atoms during 
the loss process using an appropriate sequence of laser 
pulses. Even if three resonances are experimentally ac- 
cessible the latter might yield better results as the res- 
onance with the best ratio of two-atom to single-atom 
loss can be exploited. One suitable sequence [22| which 
does not require precise control is to apply a large num- 
ber of pulses of random relative phase and approximate 
area 7r/2 each at equal intervals. Each initial state then 
spends 1/3 of the time in the resonant state, and hence 
has probability qi = exp(— ti/3rd) of failing to lose a 
pair of atoms occupying the same site where U is the 
total duration of the sequence and Td the two-atom loss 
time constant. This method will not be perfect since 
single particles are also lost from the system with some 
time constant t s and hence U cannot be chosen arbitrar- 
ily large. The probability of losing a single particle is 
pi = 1 — exp(— ti/T s ) where t s ^> r^. Both pi and qi are 
error probabilities, and their sum p\ + qi is minimized 
by choosing U = ln(r s / (3T d )) / '(1 / '3r d - 1/tA An experi- 
ment recently performed by Widera et. al |2l| measured 
Td = 1.3ms and no detectable loss of single atoms for 
resonance times up to 100ms. If we take r s = 500ms the 
above gives t\ = 19ms and error probabilities of qi = 0.8% 
and pi = 3.7%. 

In summary, the optical lattice realization of our entan- 
glement detection network contains the following stages 
at which experimental errors are likely to occur: (i) qb s 
from the implementation of Hbs] (h) Pi and qi occurring 
during the loss stage; (iii) detector errors pd in counting 
the number of singly occupied lattice sites. In addition 
(iv) the setup might lack spatial resolution in atom count- 
ing. We consider how to correct errors (i) - (iii) in Sec. IIVI 
and next study the case (iv) of no spatial resolution in 
an otherwise perfect experimental setup. 



III. ENTANGLEMENT DETECTION AND 
CHARACTERIZATION WITHOUT SPATIAL 
RESOLUTION 

We assume that the measurement of the total num- 
ber of singly/doubly occupied sites is accurate but that 
we cannot know their locations. We first show that this 
information is sufficient to detect a violation of Eq. J3J. 
Then we study how various different experimentally real- 
izable multipartite entangled states might be character- 
ized using such measurements. 
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A. Entanglement detection 

The probabilities P(j) of measuring j singly occupied 
sites in one row (2j single atoms in total) are given by 



(19) 



\A\=j B 



fc=0 v 7 I 



k\ I n — k 



where the summation indices are k — \B\,l — \Af~)B\, j = 
\A\, and A is the set of singly occupied (antisymmetric) 
sites. We form the generating function 



E^c?) - ^Eu) Tr H)}( 1 --) fc ( 1 +-) 



n—k 



1 + x 



fc=0 v 7 v 



(20) 



and let y = (1 — x)/(l + re), to obtain 



E(i - v)'(i + y) n - j P(j) = E W ! Tr H) } ■ ( 21 

7=0 fc=0 ^ 7 



J=0 

from which we find 



Tr H)} 



J=0 



^irj-Jic;.^ )(-!)' 

(22) 

Therefore, although we cannot determine the purity of a 
given subset of the row of atoms we can still determine 
average purities associated with subsets of atoms of a 
given size by measuring P(J). We will prove later (see 
Eq. (gll) that the accuracy in finding P(j) required for 

obtaining a given accuracy of Tr | p 2 k ^ | has an upper 

bound independent of n and k if no errors are present. 
The network is thus efficient in detecting the presence 
of entanglement in all pure (and some mixed) entangled 
states via the violation of Eq. Q . 

In Fig. |3| we show the probabilities P(j) and the re- 
sulting average purities for a variety of different states. 
For a classically correlated state shown in Fig. [SJi the 

values of Tr | p 2 k ^ X are monotonically decreasing with 

k showing that Eq. Q is not violated. The maximally 
entangled state shown in Fig. |3Jd has the characteristic 

that Tr {pf k) } = 1/2 for < k < n while Tr ip\ n) \ = 1 

and thus the inequalities are violated in this case. The 
cluster state shown in Fig. [3: violates the inequalities for 
all j > n/2 and therefore its entanglement is detected. 
Finally, in Fig. |3i we show a noisy cluster state which 
was affected by phase noise acting independently on each 
atom. It can clearly be seen that decoherence reduces the 
violation of the inequalities but that entanglement is de- 
tectable for small amounts of noise. 
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FIG. 3: Tr j against k (left) and P(j) against j (right) 

for n — 10. (a) Classically correlated state p = (| 0) n (0 \ n + 
| l) n (1 |")/2. (b) Maximally entangled state | (7 = 0)io). (c) 
Cluster state | {cj> = 7r)io). (d) Same as (c) with 10% dephas- 
ing decoherence. 



B. Characterization of entanglement 



The n measurable quantities Tr k = 1? • ■ • ,n 

do not provide us with enough information to determine 
an arbitrary state p. However, if it may be assumed 
that the state p is of a known form with less than n 
unknown parameters, it is often possible to determine 

these parameters from Tr 'f/ ? (fe)j- We demonstrate this 
by considering macroscopic superposition states | 7„) and 
cluster-like states | <j>j) introduced in Sec. Ill A~3l Finally, 
we will also look at product states of states of subsystems 
containing several atoms. 



1. Macroscopic superposition states 

Because the state | 7„) is totally symmetric the indi- 
vidual purities given in Eq. J7J) only depend on the size of 

the subsystem k = \B\ and thus Tr j = Tr {p%}- 

Hence, from the knowledge of Tr we can deter- 

mine the value of 7 which in principle only requires two 
of the purities. The remaining equations allow a partial 
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FIG. 4: Average purities Tr |/°( fc )| (k = 1 dotted curve, 
k — 2 dashed curve, k — 3 dash-dotted curve, k = 7 solid 
curve) for the state | n = 15 ). Since all states are pure we have 



Tr 



H)}=Tr{pU)}- 



check of the assumption that the measured state indeed 
has the form | 7„) . 



2. Cluster states 

The states | 4> n ) are parameterized by the entangling 
phase </>. The average purities as a function of 4> are de- 
picted in Fig. 2] For any value of < <fi < 2ir the states 
violate the inequalities Eq. The degree of violation 
increases with (j) until 4> = tt where the state is a clus- 
ter state and the degree of violation is a maximum. If 
(j) is increased further the state again approaches a prod- 
uct state and the degree of violation of the inequalities 
correspondingly decreases. Hence, from experimentally 

measured Tr j p^ fe ) } one can determine the phase cf> up 

to its sign. Again the over determined system (n equa- 
tions for one unknown (f) provides a check on how well 
the state fits the assumed form | </> n ). 

The effect of dephasing according to the map Eq. © 
on a cluster state is shown in Fig. |3J The average pu- 
rities decrease with increasing noise level d. Entangle- 
ment is certainly present and in principle detectable by 

our method as long as the Tr are not in descend- 

ing order, in the case shown in Fig. [S]up to d « 0.45. 
Again, the parameter d can in principle be determined 
from measuring the average purities. 



3. Products of entangled subsystem states 

Finally we give an example of a class of states where 
even though p is characterized by more than n parame- 




FIG. 5: Effect of dephasing noise on the n — 15 cluster state, 
k = 1 (dotted), k = 2 (short dashed), k = 8 (long dashed), 
k = 14 (dash-dotted), k = 15 (solid). 



ters, the associated average purities Tr on ly de- 

pend on n parameters. Consider the case where p is a 
product state of L subsystems, p = ^f^Pi, with each 
subsystem pi composed of a known number rii of atoms. 
In this case we have 



£ n(;;H4,} 



k i= k 



Since there are J2i 



n different Tr 



(23) 



in total, 



Tr | p 2 ^ | provides us with enough information to deter- 
mine the average purities Tr | p?^ fe ^ | of every subsystem. 

In particular, L = n/2 and all m = 2 is the case of only 
pairwise entanglement. Since Eq. JSJ holds for all Pi tak- 
ing rii = 2 provides a test for multi-particle as opposed 
to two-particle correlations in a given state. We note 
that unless p is pure classical multipartite correlations 
will also be detected by our network. 

We will now study how the working of the network 
is affected by experimental errors. In particular we will 
estimate how many runs are necessary to obtain the prob- 
abilities P(j) with sufficient accuracy in the presence of 
errors. 



IV. EFFECTS OF EXPERIMENTAL ERROR 

The errors introduced in Sec. Ill Cl affect the ability to 
find the purities Tr {p%} as well as the average puri- 
ties Tr associated with pi2... n - All of these er- 
rors are of one of two mathematical kinds: extra pairs of 
atoms and missing single atoms. We call these two errors 
"beam-splitter" and "detector" error respectively. Their 
respective probabilities q = qi, s + qi and p ~ Pd + Pi are 
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understood to include also errors occurring while parti- 
cles are lost from doubly occupied sites. The relationship 
Eq. fH^ between the purities of pb and the probabilities 
of an even/odd number js of singly occupied sites in 
B indicates that an experimental error, occurring with 
probability p per site, changes the result Tr {p%} by 
0(\B\p) if we do not attempt to correct for it. This ren- 
ders the measured results totally meaningless as soon as 
\B\p ~ 1, because the purity of any given state /?b is 
smaller or equal to one. However certain types of error, 
including the BS and detector errors, can be corrected 
by a suitable modification of the formulas Eqs. I|13l I22|) 

yielding Tr {p%} an d Tr This correction elimi- 



nates systematic errors, making Tr {p%} and Tr {/°(fc)l 
correct on average, but tends to amplify the random er- 
rors that are inevitable in measuring probabilities using 
a finite number of experimental runs. These random er- 
rors can in principle be made arbitrarily small for any 
p, q < 1 by increasing the number of runs, but in practice 
there is a limit because, as we will show, the number of 
runs required scales approximately exponentially in \B\p. 
We will now investigate the effects of these errors on the 
performance of the entanglement detection network both 
without and with spatial resolution. 



A. Without spatial resolution 

We assume the probabilities p and q to be the same for 
all 2n lattice sites and in the case of q for all symmetric 
atom pair states | aa) , | bb) and | ab) + \ ba) . Errors at 
different lattice sites are assumed to be uncorrelated. 



1. Beam splitter error 

Let Pexp(i) be the probability of detecting i atoms in 
an experimental run. If only BS errors are present this 
probability is given by 



Pcx P (2i) = $>Wr J,W- j {l-<lT-\ (24) 
j=0 ^ ' 

where the factor two in P oxp (2z) accounts for P(j) being 
the probability of having j antisymmetric pairs. We can 
use generating functions to invert Eq. (|24[1 



5>*P exp (2i) = £P(jV(1 -q + xq) n -\ (25) 
i=0 j=0 

leading to 



p (j) = E 



n — i\ (— q)i 



n-jj (l~q) r 



r Pex P (2i). (26) 



We now apply Eq. I|26(l to a subsystem B and substitute 
this into Eq. giving 

Tr {pD = E (l^V ' B (-l) ifl Pex P (2iB), (27) 



■i B =0 



<1 



where is refers to the number of atoms detected in B. 
This expression is then averaged over all B of size \B\ = k 
to give 



Tr {p(k)} =E Afe - p -p( 2z )> 



(28) 



i=0 



with 



Ai 



n -1 fc 



1-9 



k-l 



IJ \ k — I 



(29) 

Hence, using Eq. I|29|l instead of Eq. (|22|) corrects all 
the systematic error caused by an imperfect BS. We are 
still left with the inherent random error associated with 
the measurement of the probabilities P oxp (2i), which is 
reduced by increasing the number of experimental runs. 

Because of this random error the estimate of Tr |p( fe ) j 

obtained from N experimental runs (each using one pair 
p® p) has the correct mean but a nonzero standard devi- 
ation y^Vk/N, where this defines Vk{p, q, p)\ hence 0(Vk) 
runs are necessary for meaningful results. Note that in 
general Vk > even when p = q = 0, as it includes the 
inherent quantum uncertainty as well as that added by 
experimental error. In the case of BS error 



V k = E p ex P (i)4i- (Tr {pf k) }) <max i (A| 



< 



1+q 
1-q 



2/,- 



(30) 



where the approximation is valid for > 1, g < 1. The 
bound Eq. I|3U|) proves that the number of runs required 

to obtain meaningful estimates of Tr | p^.) } is reasonable 
for k < 1/g, however large n is. 

We numerically computed the worst case by maximiz- 
ing Vk with respect to Tr |p(fe)| subject only to Eq. (0 

and compare it to cluster states | (4> — 7r)„) in Fig. 
The results confirm the analytically found exponential in- 
crease of Vk with q in the worst case. For the cluster state 
Vk increases only slowly with q for small k < 1/q while 
for k > 1/q we find approximately exponential growth of 
Vk with q. We also computed the variances for maximally 
entangled states and found that they are quite close to 
the worst case shown in Fig. [7^,. Therefore one may in 
an experiment generally not expect the variances Vk to 
be much smaller than our worst case results. Thus only 
BS errors up to q = l/k are acceptable and yield reliable 
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FIG. 6: Variance Vk against BS error, for k = 1 (dotted), 
fc = 4 (short dashed), k — 7 (long dashed), fc = 11 (dash- 
dotted), ft = 15 (solid) curve, and n — 15. (a) shows the 
result for the worst case and (b) for a cluster state. 




FIG. 7: Worst case variance against n with BS error q — 
1/k, for k — 2 (diamonds), A; = 4 (triangles), k — 7 (hollow 
squares), k = n (solid squares). The dotted horizontal line 
is the analytic bound Vk < exp(4fcq). The solid curves are 
drawn to guide the eye. 



results in a reasonable number of runs for all < k < n. 
However, as shown in Fig. for q = 1/k, the average 
purities with k <C n will be determined much more accu- 
rately than those with k w n and should thus be prefer- 
entially used for determining parameters characterizing 
the measured state. 

We finally note that if J is fluctuating from run to run 
and qi, s hence becomes a random variable error correction 
is still possible. In this case we have to replace Eq. H24|) 
by 



P exp (2i) = J2 P U) 



n-j 



dJf(J)q(jy-i(l-q(J)T-\ (31) 



where /(J) is the probability density function of J and 
q (J) denotes the BS error as a function of J. The result- 
ing system of linear equations Eq. (|31[) can be treated 



using the methods introduced above. 

2. Detector error 

We now assume that only detector errors are present 
and each atom has a probability p of failing to be de- 
tected. In this case P exp («) is related to P(j) via 

'2/ 



('■) 



E p u) 

j=i/2 



2j- 



(32) 



We can again solve this equation by methods similar to 
those used in Sec. IIV A II (though this time it is not a 
true inverse because the system is over-determined |23j) 
to obtain 



P(j) 



E 

i=2j 



(-P) 



i-2 j 



2jJ (1-pY 



"-^cxp(^) ) 



(33) 



where non-integer values of j are discarded. By combin- 
ing Eqs. H22II33I) we obtain Tr m terms of P CX p(i)- 

For the remaining random error measured by Vk we ob- 
tain the upper bound 



Vk < 



l+p 
1-p 



a 8np 



(34) 



Numerical calculations for n < 15 confirm this exponen- 
tial growth of Vk with p at np ~ 1. The results are 
shown in Fig.|Sj This time, however, Vk is typically much 
smaller than the analytic bound, e.g. for the n = 15 clus- 
ter state, fitting Vk oc exp(/3np) gives (3 ~ 2. However, 
the exponential growth with n implies a practical limit 
of n ~ 1/p for any k contrary to the case of BS errors. 
This scaling can be improved by using the least squares 
method to handle the over-determined linear system of 
equations Eq. i',i'2\) . We do not have an analytic bound 
analogous to Eq. H34|) for the least squares method but 
numerical calculations for the cluster state and the worst 
case result in significantly lower values for Vk than those 
obtained from Eq. I|33l) . Most importantly as shown in 
Fig. it appears that the scaling of T4 becomes expo- 
nential in kp rather than np similarly to the case of BS 
errors. This implies that an error of p < 1/k is accept- 
able for obtaining meaningful estimates of Tr | p 2 ^ | in 
a reasonable number of experimental runs. 

In an actual experiment both BS and detector errors 
will be present and first correcting for the detector error 
p using Eq. I|33() . then substituting the resulting P(j) for 
-Fexp(2j) in Eq. H29|) to correct for the BS error q yields 
a combined analytical error bound of 



V k < 



l+p 



1 



P 



1 + q 
1-q 



2 k 



a Snp-\-4:kq 



(35) 



According to our numerical results using the least squares 
method this bound can be improved requiring only p,q < 
1/k for obtaining sufficiently small errors. 
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FIG. 8: Variance against detection error without spatial reso- 
lution using Eq. (I33H . for k — 1 (dotted), k — 4 (short dashed), 
k — 7 (long dashed), k = 11 (dash-dotted), k = 15 (solid), 
and n = 15. (a) shows the result for the worst case and (b) 
for a cluster state. 
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FIG. 9: Worst case variance against n with detection error 
p — l/(2k) using least squares, for k — 1 (diamonds), k = 2 
(triangles), k — 3 (hollow squares), n (solid squares). The 
solid curves are drawn to guide the eye. 

B. With spatial resolution 

If spatial resolution is available then Tr {p%} , not just 

its average Tr over all subsets B of size k, be- 

comes accessible to measurement. This allows us to do 
some state characterizations which would otherwise be 
impossible and also introduces extra redundancy. As we 
will show below this does not affect the tolerance to BS 
errors but we will find that the detector error tolerance 
improves to p 2 ~ 1/| B\. Imperfections in the spatial 
resolution, however, will lead to additional errors in de- 
termining Tr {ps}- 

1. Beam splitter error 

The variance of Tr {p%} due to BS errors can be di- 
rectly inferred from Eq. (|27[l . where only the atoms in B 



are counted. Using the same methods as in Sec. IIV A II 
we find an upper bound for Vb (defined analogously to 
V k ) given by V B < ((l+g)/(l-<?)) 2|B| « exp(4|B|g). The 
BS error tolerance (for fixed k = \B\) is thus independent 
of whether spatial resolution is available or not. 

2. Detector error 

The situation is different for the detector error. Each 
antisymmetric pair contains two atoms only one of which 
needs to be detected to know that it was antisymmetric. 
Therefore the effective error probability becomes p 2 . The 
resulting formula is 

Tr{p|}=EMr(^)Vf xpW , (36) 

where P^ p (i) is the probability of measuring i antisym- 
metric sites in subset B. The variance bound is given by 
V B < ((l+p 2 )/(l-p 2 )) 2 ^ « exp(4|i?|p 2 ). 

3. Imperfect spatial resolution 

There is a new type of error to consider as the spatial 
resolution itself will not in practice be perfect. If we 
let f(x,y) be the probability of finding at position x a 
particle which is actually at position y, then 

1 k 

P C * P (A) =EE^)I1 f(M+k) - Bi)/(^ W , S<)P(B), 

i£B ? ^ ' i=0 

(37) 

where A is the "set" (its elements are not necessarily 
distinct) of observed atom positions, and P exp (A) is the 
experimental probability of observing atoms exactly at 
these positions A. The set B denotes the antisymmetric 
sites and the probability of having antisymmetric sites at 
positions B is P(B), k = \B\ = \A\/2, and ? runs over 
all (2k)\ permutations of the 2k atoms in A. The i-th 
element of B and A are written as Bi and Ai, respec- 
tively. The symmetry factor s(A) stands for the num- 
ber of permutations ? which leave the ordered lists of 
atoms {Ai} invariant (e.g., s({l, 1}) = 2, s({l,2}) = 1), 
and is needed because our summation runs over differ- 
ent ordered lists {Ai} of the same set A. This is an 
over-determined linear system, and just as in the case of 
detector error, we can either explicitly solve it by dis- 
carding some of the equations or apply the least squares 
method. As an explicit solution we can e.g. use 

P (B) = E 4^ mr\B u Ai)\ Pexp(A), (38) 

{Ai} \i=l / 

where we define Bk+i = Bi. The sum runs over all or- 
dered lists {Ai} and / _1 is the matrix inverse of /, that 
is f^ 1 (y,x)f(x,z) = 5 XZ . 
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FIG. 10: Variance of Tr {p|} as a function of the standard 
deviation a, for B = {2} (dotted), B = {2,3} (dashed), B = 
{1,3} (dash-dotted), B = {1,2,3} (solid), and n = 4. (a) 
and (b) are worst case variances while (c) and (d) are for a 
maximally entangled state, (a) and (c) each show four almost 
coincident curves obtained from Eq. 1381 while (b) and (d) are 
found by the least squares method. 



macroscopic superposition states. We investigated the 
experimental realization of the network in an optical lat- 
tice and identified lack of spatial resolution, errors in the 
BS operation and imperfect atom detection as the main 
sources of error. We showed that even in the absence 
of spatial resolution entanglement can be detected. In 
cases where the entangled state is characterized by a few 
parameters only we found that these can often be deter- 
mined from the measurement results. We also studied 
the influence of BS errors occurring with probability q 
and detection errors with probability p and concluded 
that with small numbers of experimental runs entangle- 
ment can be detected between k atoms as long as kp < 1 
and kq < 1. Finally we showed that for obtaining pu- 
rities Tr {p%} of subsets B rather than average purities 
with reasonable experimental effort a spatial resolution 
of a w A is necessary. 

The results obtained in this work show that unam- 
biguous multipartite entanglement detection in optical 
lattices is possible with current technology. This has 
not yet been achieved experimentally 6]. Furthermore 
it will even be possible to determine some of the charac- 
teristics of entangled states created in these experiments 
without the requirement of performing spatially resolved 
measurements. Our network is thus a viable alternative 
to detecting entanglement via witnesses or full quantum 
state tomography. 



As before by performing numerical calculations using 
the least squares method we find much lower variances 
than with Eq. I|38|) . An example is shown in Fig. ITUI 
where we plot the variance of Tr {p%} due to a Gaussian 
position error of the form 

x+\/4 

f(x,y) = -jL- J dzexp[-(z-y) 2 /(2a 2 )} (39) 

x-X/4 

where a is the standard deviation and A is the wave 
length of the laser creating the optical lattice. The cor- 
responding lattice spacing is A/2. The results obtained 
from Eq. (|38|) shown in Fig. HUfa . c require a resolution 
of a < A/2 whereas the least squares method shown in 
Fig.llUb. d yields reasonable variances Vb for spatial res- 
olutions up to a < 3A/2. However, due to the exponen- 
tially large number of possibilities for A the least squares 
calculation becomes intractable for large n. 

V. CONCLUSION 

We discussed the detection and characterization of 
multipartite entanglement in optical lattices with the 
quantum network introduced in under different ex- 
perimental conditions. We first described how the net- 
work can be implemented in ideal experimental condi- 
tions and showed that it allows to characterize a number 
of important classes of states like cluster-like states and 
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APPENDIX A: THE BEAM SPLITTER 
OPERATION 

Since the BS only couples two lattice sites in rows I and 
II of each column (see Fig. |2J we consider a single such 
pair and omit the column superscript j in this section. 
For U = 0, we obtain from the Heisenberg equations for 
the operators a = a,b 

aj{t) = cos( Jt)ai — i sin( Jt)au, 
otn(t) = cos( Jt)otn — isin(Jt)ai. (Al) 

Hence, applying Hbs f° r a time tb s = ?r/(4J) implements 
a perfect pairwise BS. 

Initially the atom pair is in a state of the form p® p, 
where p is a single qubit state and hence has a spectral 
decomposition of the form p = Ai | c) (c I+A2 \d) (d\, with 
Ai + A 2 = 1, I c) = c t I vac), and \d) = S \ vac). Here 
c\ d) are linear superpositions of a' , w with coefficients 



12 



depending on p. Therefore we can write 

p <8> p = \\ | cjc/j) (c/c// I + Ag j did n ) (didu I (A2) 
A1A2 



+ 



2 

A1A2 



(I c/d//) + I djcu))({cidii I + (rf/c// |) 
(I c/d//) - I d I c H )){(c I d II I - (d/c// |), 



which is a classical mixture of a symmetric state with 
probability P + = 1 — A1A2 and an antisymmetric state 
with probability P_ — A1A2. After the BS the resulting 
state p' = exjp(iHBst)p <8> pexp(—iHBst) is given by 

p' = A 2 |$i)($i|+A 2 |$ 2 )($ 2 | (A3) 

+ A1A2 I $3) (*3 I + A1A2 I $4) ($4 I , 

where | $1) = {c\c\ + c\ lC \j) | vac) /2, | $ 2 ) = {d\d\ + 

4/4/) I vac > A I *s) = (44 + 4/4x) I vac ) are 
states with double occupancy in one row and an 

empty site in the other row while | $4) = (c\d\j — 



in each row. Hence, after the BS we will find a doubly 
occupied site with probability 1 — A1A2 = P+ while two 
singly occupied sites result with probability A1A2 = P-. 

However, in practice it is impossible to completely turn 
off the interaction U which may result in a symmetric 
state failing to bunch. We consider any symmetric state 
I VPg) and because Hbs + Hjj acts only on the row indices, 
not the internal states, and is symmetric between the two 
rows the probability of failure to bunch is given by 



qbs = 



(m s \e l(HBS+Hu)t \y s 
16J 2 



(A4) 



+ 



u 2 



16J 2 + £/ 2 ' 



The optimal choice for the BS time is = 
7T/V16J 2 + U 2 for which q bs = C/ 2 /(16J 2 + U 2 ). If 
the hopping term is not controlled perfectly accurately 
but fluctuates by SJ around a mean J we set tb s — 



4) I vac) /^/2 is a state with a singly occupied site ir/^WJ 2 + U 2 and obtain Eq. (|18fl . 
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